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ABSTRACT 

Aims. We investigate the binary star T Tauri South, presenting the orbital parameters of the two components and their individual 
masses. 

Methods. We combined astrometric positions from the literature with previously unpublished VLT observations. Model fits yield the 
orbital elements of T Tau Sa and Sb. We use T Tau N as an astrometric reference to derive an estimate for the mass ratio of Sa and Sb. 
Results. Although most of the orbital parameters are not well constrained, it is unlikely that T Tau Sb is on a highly elliptical orbit 
or escaping from the system. The total mass of T Tau S is rather well constrained to 3.0*2H M G . The mass ratio Sb:Sa is about 0.4, 
corresponding to individual masses of Ms a = 2.1 ± 0.2 M Q and Msb = 0.8 ± 0.1 M Q . This confirms that the infrared companion in the 
T Tauri system is a pair of young stars obscured by circumstellar material. 
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1. Introduction 

T Tauri was discovered in 1852 by Hind during observati ons of 
the br ight variable nebula N GC 1555 (H ind's nebula) dBertoull 
1984). Almost a century later. lJoyl d 19451) defined a new class of 
variables, which contained eleven stars initially. The members 
of the class were later identified as low-mass pre-main-sequence 
objects. Joy named these sources 'T Tauri variables', because 
T Tau not only shows the typical features, but also is among the 
brightest and best known objects of this group. 

However, as observational techniques improved, the charac- 
ter of T Tau turned out to be more peculiar th an prototypical. On e 
of the most surprising findings was made bv lDvck et alj £l982). 
They observed T Tau using one-dimensional speckle interferom- 
etry and found an infrared companion with a separation of less 
than one arcsecond. However, neither the relative position of the 
two sources nor the identification of the vis ually bright compo- 
nent could be derived. Radio observations dCohen et alj fl982) 
and subsequent astrometric comparisons between the radio and 
the optical data revealed that the infrared companion was lo- 
cated 077 ± 0".' 1 dde Vest! 19821) or 076 + 07 1 dHanson etalj 19831) 
south of the visually bright component. The southern component 
dominates the system longward of the near-infared wavelength 
regime, but no optical c ounterpart could be id entified down to 
19.6 mag in the V-band dStapelfeldt et alj|1998l) . 

Explanations of the physical nature of the enigmatic in- 
frared companion ranged from a true pr otostar still embedded in 
its contracting envelope dBertoutlll983l) . to a FUOr- like source 
with the observed rapid brightening due to accretion dGhez et alj 
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119911) . to a normal T Tauri star coeval with its primary. In the lat- 
ter case, the star has to be obscured either due to a s pecial view- 
ing g eometry through the ci rcumstellar material dCalvet et al.1 
1994) or enhanced accretion dKoresko et al.ll 1997b . 

Considerable observational and theoretical effort has been 
expanded to understand this source. Perhaps the most important 
ingredient in solving this puzzle came surprisingly: the resolu- 
tion of the s outhern compa nion into two sources using speckle 
holography (Koresko 2000). With a separation of ~ 0705, the 
two components were named T Tau Sa and T Tau Sb. T Tau Sb 
has been identified as a heavily extincted, act ively accreting 
pre-m ain-sequence star with a spectral type Ml (Duchene et al. 
|2002|) . while the spectrum of T Tau Sa is featureless except for a 
strong Bry emission line, leaving the nature of the object uncer- 
tain. 

The northern o f the two radio sources found by 
ISchwartz et all d 19841) is easy to identify with T Tau N, but it 
was not clear which of the two componen ts of T Tau S should be 
ident ified with the s outhern radio source. Uphnston et al. (2003) 
and lLoinard et al.1 d2003l 120051) associated the radio source 
with T Tau Sb. This would require that the orbital motion of 
T Tau Sb around T Tau Sa shows a dramatic change between 
1995 and 1998, probab ly caused by th e ejection of Sb from 
the T T au S system dLoinard et al.ll2003h . However. iFurlan et al.l 
(2003) found that T Tau Sb and the radio source have distinct 
paths, and suggested a fourth object as the source of the south- 
ern radio emission. After further refinement of the picture, it has 
been concluded that the radio source may be connected wit h, 
but not identical to t he infrared source dJohnston et al.ll2004allbt 
lLoinard et ai1l2007al) . 

For completeness, we note that the detection of an- 
other source 0727 no rth of T Tau N has been reported by 
iNisenson et al.l d 1 985b . They identified a second component in 
the visual wavelength regime, which was later detected in the 
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K-band by Maihara & Kataza ( 1991). Howeve r, this source was 
not de tected by iGorham et"ai] (1 1992b nor by IStapelfeldt et al.l 
(119981) . If the additional source exists, probable explanations for 
the non-detections are that the source is variable either intrin- 
sically or due to varying extinction. The latter i s not unlikely 
considering the complex environment of T Tau dSolf & Bohml 
[19991: iHerbst etai]2007h . 

The presence of the companion Sb at such a small separa- 
tion offers the possibility to trace the orbital motion of T Tau Sa 
and Sb around each other in a reasonable time and with no in- 
fluence of T Tau N in first order approximation. This can lead 
to the determination of the dynamical mass of T Tau S and - 
with the help of the common orbit of the southern pair around 
the northern component - the individual masses of T Tau Sa and 
Sb. Such information is essential to judge the different models. 
In addition, evolutionary models of pre-main-sequence objects 
may benefit from the direct measurement of stellar masses, since 
only a few young stellar objects have masses known independent 
of th eoretical assumptio ns. 

iDuchene et al.l d2006) made such a dynamical analysis of the 
T Tau S system using measurements in both the near-infrared 
and radio. A simultaneous fit of the Keplerian orbit of T Tau Sb 
around T Tau Sa and the orbital motion of T Tau Sb around 
T Tau N lead to a dynamical mass of 3.34+0.36 M for the south- 
ern binary and a mass ratio Sb:Sa of 0.18 ± 0.04. This result is 
consistent with earlier determinations of the individual masses 
dJohnston et al. 2004a b). The improved accuracy allows us to 
conclude that T Tau Sa with 2.73 + 0.3 1 M Q may be a Herbig Ae 
star wh ose light is extincted by a highly inclined circumstel- 
lar di sc dKasper et alJl2002tlDuchene et alJl2005l;lDuchene et alj 
120061) . The derived mass of T Tau Sb is 0. 61 +0.17 M , in good 
agreement with its early M-type spectrum. lDuchene et alj d2006) 
modeled the orbital motion of T Tau S around N by polynomi- 
als in x and y. This does not allow a determination of the sys- 
tem mass, and prevents the derivation of the dynamical mass of 
T Tau N. 

In this paper, we further refine our knowledge of the or- 
bits and thus of the individual masses. We restrict ourselves 
to measurements made in the near-infrared, since the iden- 
tification of the radi o sources with T Tau Sb is uncertain 
dLoinard et al.ll2007bb . A further complication comes from the 
fact that T Tau Sa has not been detected in radio observations. 
The position of Sb relative to Sa can only be computed indirectly, 
using the position of T Tau N and the estimated parameters of the 
orbits. 

We also explore the confidence limits for the resulting orbital 
elements, a crucial step when interpreting such highly degener- 
ate orbits as those in the T Tau triple system. 

The data from the literature, as well as the reduction and 
analysis of our own measurements, are presented in Section [2] 
In Section[3] the fit of the orbital elements and the determination 
of the mass ratio is explained. The reliability of the fit and its 
implications for our current knowledge of the T Tau system is 
discussed in Section|4] In Section[5]we summarize our results. 



2. Observations 

2.1. Literature 

T Tauri has been observed in the near-infrared with various tele- 
scopes and instruments since the discovery of its binary nature. 
Most of the measurements used in this work thus were taken 
from the literature (see Tables Q] and |2). The tables include all 
statistical and systematic errors reported by the authors in the pa- 



Table 1. Astrometric Measurements of T Tau N 
solved) 



S (S unre- 



Date (UT) 



d [mas] 



p.a. n 



Source 



1989 Dec 10 

1990 Nov 10 

1991 Nov 19 

1992 Feb 20 

1992 Oct 11 

1993 Nov 25 

1993 Dec 26 

1994 Sep 22 
1994 Oct 19 
1994 Dec 25 

1997 Nov 17 

1998 Nov 2 

1999 Nov 23 
1997 Dec 4 
1997 Dec 6 

2000 Feb 20 



720 ± 14 
716 + 14 
702 + 14 
710+14 
692 ± 15 
701 ± 14 

690 ± 14 
689 ± 14 
689 ± 14 
718 ± 9 
709 ± 6 
704 ± 6 
699 ± 6 
685 ± 13 
698 ± 13 

691 ± 2 



175.5 : 

176.2: 
176.4: 
176.7: 
177.0: 
177.0: 

176.0 : 
175.9: 
178.2: 
177.6: 
180.2: 
180.7 : 
181.0: 
179.5 : 

179.1 : 
181.7: 



1.0 
1.0 
1.0 
1.0 
1.4 
1.0 

1.1 

1.2 

1.1 

0.7 
0.5 
0.5 
0.5 
1.0 
1.0 
0.8 



Ghez et al 
Ghe z et al 
Ghez et al 
Ghez et al 
Ghez et al 
Ghez et al 
Ghez et al 
Ghe z et al 
Ghez et al, 
Roddi er et al 
Roddier et al 



(1995) 
(1995) 
(1995) 
(1995) 
Q995) 
(1995) 
(1995) 
(1995) 
(1995) 
2000) 
2000) 



Roddier et al. (2000) 
Roddier et al. (2000) 
White & Ghez (2001) 
Whit e & Ghez (2001) 
This work 



pers or in follow-up studies. The data of the observations that did 
not resolve the southern binary are listed in TableQ] Table[2]hsts 
observations that did resolve T Tau S and gives the p ositions of 
compo nent Sa relative to N and Sb, respectively. Duchen e et al.l 
(2006) give only positions of Sb relative to N and Sa; the po- 
sition of Sa relative to N has been computed by us, taking the 
propagation of errors into account. 



2.2. ALFA 



iKohler et al.l d2000al) presented results of observations obtained 
in February 2000 with the adaptive optics system ALFA at the 
3.5 m telescope on Calar Alto, Spain. The H photometric band 
was used to take advantage of the smaller diffraction limit com- 
pared to the K band. T Tau N served as PSF-reference to de- 
convolve T Tau S using speckle interferometric techniques. This 
allowed t hem to resolve the Sa-Sb binary, and confirm the dis- 
covery of lKoreskol (120001) . 

We revisited these observational data, in order to measure 
the position of component Sa relative to N. In our new analysis, 
we used observations of the star PPM 1 19798, taken shortly af- 
ter the images of T Tauri, to deconvolve the triple system. We 
again applied our own speckle pro gram; for a detailed descrip- 
tion see, e.g.. IKohler et al.l d2000bl) . The resulting visibility ap- 
pears in Fig.Q] It clearly shows the characteristic fringe pattern 
of the N-S binary. However, the signal of the Sa-Sb pair is ob- 
scured by random fluctuations. One has to keep in mind that 
T Tau N contributes more than 90 % of the flux of the system 
in the H band, resulting in correspondingly high speckle noise. 
Therefore, it was not possible to measure the relative position of 
T Tau N and Sa, but only the position of the combined light of 
Sa and Sb relative to N. The latter is listed in Table [1] We then 
computed the offset of T Tau Sa from the center of light of Sa 
and Sb, using the known position of Sa relative to Sb an d the 
flux ratio of Sb/Sa of 0.259 + 0.011 (iKohler et al.l l2000a). The 
results are given in Table [2] 

2.3. N AOS/CONIC A 

We also report on observations with NAOS/CONICA (NACO 
for short), the adaptive optics, near-infrared came ra at the ESQ 
Very Large Telescope on Cerro Paranal, Chile dRousset et al.l 
2003; iLenzen et all 120031) . They are marked "this work" in 
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Date(UT) T Tau N - Sa TTauSa-Sb Source 

d [mas] P. A. [°] d [mas] P. A. [°] 



1997 Oct 12 


693 ± 


12 


178.6 ± 


0.9 


51 


± 


9 




218 


+ 


8 


Duchene etal. (2006) 


1997 Dec 15 


— 




_ 




53 


+ 


9 




225 


+ 


8 


Koresko (2000) 


2000 Feb 20 


686.3 ± 


2.4 


180.5 ± 


0.8 


79 


+ 


2 




253 


+ 


2 


This work / Kohler et al. (2000a) 


2000 Nov 19 


702 ± 


5 


180.8 ± 


0.2 a 


92 


+ 


3 




268.1 


+ 


1.6 a 


Duchene et al. (2002) 


2002 Oct 30 


697.6 ± 


2.1 


183.3 ± 


0.3 


106.5 


+ 


2. 


.6 


284.5 


+ 


1.4 


Schaefer et al. (_2006) b 


2002 Nov 20 


693 ± 


4 


183 ± 


1 


94.3 


+ 


3. 


.3 


278 


+ 


1 


Mavama et al. (2006) 


2002 Dec 13 


695 ± 


7 


183.3 ± 


0.7 


108 


+ 


1 




284.9 


+ 


0.9 


Duchene et al. (2005) 


2002 Dec 15 


693.5 ± 


1.4 


182.3 ± 


0.3 


106.4 


+ 


1. 


.6 


283.9 


+ 


0.4 


This work 


2002 Dec 20 


— 




_ 




110 


+ 


4 




282 


+ 


2.1 


Beck et al. (2004) 


2002 Dec 24 


691 ± 


12 


181.7 ± 


1.1 


108 


+ 


5 




287.8 


+ 


2.6 


Furlan et al. (2003) 


2003 Nov 18 


700.4 ± 


2.9 


183.2 ± 


0.6 


116.5 


+ 


3. 


.1 


291 


+ 


1.6 


Schaefer et al. (2006) 


2003 Dec 12 


698 ± 


5 


181.9 ± 


1.2 


118 


+ 


2 




288.6 


+ 


1.1 


Duchene et al. (2005) 


2003 Dec 12 


696.2 ± 


2.0 


182.8 ± 


0.1 


113.3 


+ 


1. 


.0 


288.5 


+ 


0.8 


This work 


2004 Nov 23 


690 ± 


3 


186 ± 


1 


100 


+ 


2 




298 


+ 


1 


Mavama et al. (2006) 


2004 Dec 9 


693.2 ± 


1.2 


183.5 + 


0.1 


118.9 


+ 


1. 


.1 


295.9 


+ 


0.6 


This work 


2004 Dec 19 


689 ± 


3 


183.6 ± 


0.4 


116 


± 


4 




294.1 


+ 


1.4 


Duchene et al. (2006) 


2004 Dec 24 


698.5 ± 


2.2 


184.7 ± 


0.3 


119.6 


+ 


2. 


.8 


297.2 


+ 


1.4 


Schaefer et al. (2006) 


2005 Feb 9 


704.0 ± 


1.6 


184.4 ± 


0.1 


120.7 


+ 


1. 


.9 


300.4 


+ 


0.9 


Schaefer et al. (2006) 


2005 Mar 9 


702.8 ± 


2.4 


184.4 ± 


0.2 


120.3 


+ 


2. 


.9 


300.6 


+ 


1.4 


Schaefer et al. (2006) 


2005 Mar 24 


703.7 ± 


0.9 


184.3 ± 


0.1 


121.5 


± 


1. 


.2 


300.4 


+ 


0.6 


Schaefer et al. (2006) 


2005 Oct 20 


703.3 ± 


1.2 


184.9 ± 


0.1 c 


123.7 


± 


1. 


.3 


303.3 


+ 


0.6 


Schaefer et al. (2006) 


2005 Nov 13 


694 ± 


2 


184.4 ± 


0.2 


119 


+ 


1 




300.6 


+ 


1.0 


Duchene et al. (2006) 


2005 Dec 6 


703.0 ± 


5.7 


185.8 ± 


0.5 


121.9 


+ 


6. 


.8 


304.8 


+ 


3.2 


Schaefer et al. (2006) 


2006 Oct 1 1 


694.4 ± 


0.9 


184.9 ± 


1.0 


125.1 


± 


0. 


.7 


304.7 


+ 


1.0 


This work 


2007 Sep 16 


693.6 ± 


1.8 


185.0 ± 


0.2 


127.0 


+ 


0. 


.5 


308.2 


+ 


0.3 


This work 



a Corrected value according to lDuchene et alj {20061) 

b This is a re-analysis of the data presented in Beck et al. ( 2004) using a more sophisticated PSF fitting routine. 

c 90° subtracted from the position angle in Schae fer et al.1 feOOq) to correct for the orientation of the camera (Schaefer, priv. comm.). 




Fig. 1. Result of the speckle-observation of T Tauri on February 
20, 2000. Shown are modulus (left) and phase (right) of the com- 
plex visibility (the Fourier-transform of the image). The N-S bi- 
nary pair causes the fringe pattern in the modulus and the steps in 
the phase. The fringe pattern due to the binary nature of T Tau S 
is too weak to be recognizable. 

Table[2] The NACO observations in 2002 and 2003 were carried 
out in the course of the NACO Guaranteed Time Observations 
(PI Tom Herbsfl), while the observations in 2004, 2006, and 
2007 were regular proposals for open time. We use only imag- 
ing observations in the K or K s photometric band for the orbit 
determination. 

The NACO images were sky subtracted with a median sky 
image, and bad pixels were replaced by the median of the clos- 



1 iHerbst et alJ d2007h give astrometric data derived from the obser- 
vations in 2002. However, their analysis was carried out independently 
from ours. 



est good neighbors. Finally, the images were visually inspected 
for any artifacts or residuals. Figure [2] shows an example of the 
results. 

The Starfinder program dDiolaiti et alj [2000) was used to 
measure the positions of the stars. The positions in several im- 
ages taken during one observation were averaged, and their stan- 
dard deviation used to estimate the errors. To derive the exact 
pixel scale and orientation of the detector, we took images of 
fields in the Orion Trapezium during each observing campaign. 
The measured posit ions of the stars were compared with the co- 
ordinates given in McCaughrean & Stauffer] d 19941) by the as- 
trometric software 

ASTROM3- 

The calibrated separations and 
position angles appear in Table [2] 



3. Orbit Determination 

3.1. The Orbit of T Tauri Sa/Sb 

We estimated the orbital parameters of the Sa-Sb pair by fitting 
orbit models to a ll observations l isted i n Table [2] except the po- 
sitions given by Mava ma" et al.l (|2006). These two data points 
show significant offsets from all the other measurements (Fig. [6] 
and|7]l. Since we do not know the reason for this discrepancy and 
have no way to correct it, we excluded these points from the fit. 
To avoid any systematic errors due to uncertainties of the dis- 
tance, the separations are expressed in units of milli-arcseconds 
within the fitting program. Only the final result was converted to 
AU. 

To find the orbital fit with the mi nimum y 2 , we employed 
a method similar to that outlined in Hilditch (2001): a grid- 



see http://www.starlink.rl.ac.uk/star/docs/sun5.htx/sun5.html 
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Fig. 2. Image of the T Tau system taken with NACO in the K s 
photometric band on 2006 Oct 11, shown with a logarithmic 
scale. 



search in eccentricity e, period P, and time of periastron Tq- 
At each grid point, the Thiele-Innes elements were determined 
by a linear fit to the observational data using Singular Value 
Decomposition. From the Thiele-Innes elements, the semima- 
jor axis a, the angle between node and periastron u>, the position 
angle of the line of nodes £2, and the inclination i were computed. 

We decided to scan a wide range of parameter values: 200 
points within < e < 2, 250 points within lOyr < P < 3100yr, 
and initially 100 points for To distributed over one orbital pe- 
riod. After the initial scan over To, the best estimate for To was 
improved by re-scanning a narrower range in To centered on the 
minimum found in the coarser scan. The range was reduced by 
a factor of 10, but the number of grid points was held constant. 
This grid refinement was repeated until the step size was less 
than one day. 

For every point in the e-f-plane, the set of orbital elements 
resulting in the minimum y 2 was stored, to allow later analysis 
of x 2 as a function of e and P. The parameter giving the time 
of periastron Tq is periodic: we obtain the same orbit if we add 
one orbital period to To. Furthermore, a change in To only shifts 
the position at a given time within the orbit. A value of To that 
has a small offset from the optimal To will therefore shift the 
expected positions away from the measured positions and result 
in a larger^- 2 . For these reasons, we are confident that our grid 
search found the globally optimum To and did not store all the 
sets of orbital elements for non-optimum To. 

The resulting ^-minimum is very broad (Figs. |3]— 0- For 
example, Fig. [3] shows that orbits with periods between 30 and 
more than 1000 yea rs are within Ay 2 < 1 of the best orbit (the 
lcr confidence limit, iPress et aljfl992h . The situation for the ec- 
centricity e is similar (Fig. [4]): Eccentricities between about 0.4 
and 0.8, and even some unbound orbits with e > 1 are within 
A* 2 < 1. 

We improved the results of the grid-s earch with a Leve nberg- 
Marquardt x 2 minimization algorithm dPress et al.l [1992) that 



154 ; 

153 r 
152 r 
151 r 
150 r 




149 E , , , , , : 

10 100 1000 

Period [years] 

Fig. 3. The result of the first step in our orbit-fitting procedure. 
Shown is the minimum x 2 for a given period, i.e. we hold the 
period fixed and search for the optimum set of the remain- 
ing 6 orbital elements. Solutions with unphysically high system 
masses ( > 250 Mq) have been excluded. The ripples at periods 
^ 200 years are artefacts caused by the grid search. Note that the 
X 2 shown here has not been divided by the number of degrees of 
freedom. 



154 ; 

153 r 

152 '- 
l x : 
151 r 

150 r 




149 E , , , , i , , , , i , , , , i , , , ^ 

0.0 0.5 1.0 1.5 

Eccentricity 

Fig. 4. Like Fig. [3] but showing the minimum^ 2 as function of 
eccentricity instead of period. Note that some orbits with e > 1 
(i.e. unbound orbits) have a x 2 not much higher than the mini- 
mum. 



fits for all 7 parameters simultaneously. The simple-minded ap- 
proach would be to use the orbital elements with the minimum 
X 2 found with the grid-search. However, given the structure of 
the x 2 -plane (Fig. |5), there is always the danger that the al- 
gorithm converges on a local instead of the global minimum. 
To avoid this, we decided to use all orbits resulting from the 
grid-search as starting points and carried out 250 x 200 runs of 
the Levenberg-Marquardt algorithm (the grid spans 250 periods 
and 200 eccentricities). The orbit with the globally minimum^ 2 
found in this way is shown in Figs. [6] and [7] and its elements are 
listed in Table[3] 
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Table 3. Parameters of the best orbital solution for Sa - Sb. 



2.0 
1.5 - 

£■ 
o 

o 
o 
LU 

0.5 - 

0.0 _ 
10 




100 

Period [years] 



1000 



Fig. 5. Contour plot of the minimum^ 2 as function of period and 
eccentricity. The contour^ 2 = 151 marks the 68% confidence 
intervals for the period or eccentricity (cf. Figs.[3]and|4]l. 



Orbital Element 



Date of periastron To 

Period P (years) 
Semi-major axis a (mas) 
Semi-major axis a (AU) 
Eccentricity e 

Argument of periastron oi (°) 
P. A. of ascending node Q (°) 
Inclination (' (°) 
System mass M s (M Q ) 
reduced x 2 



Mass ratio M sb /M Sa 
Mass of Sa M Sa (M G ) 
Mass of Sb M sb (M Q ) 



Value 



2451091 
(1998 Oct 4) 

Q-J +115 
-55 

201 +430 
zul -no 

30 +63 
57 +a2 ° 

J ' -0.22 
+12.0 
-30.0 

283.2 
55 1 +2 -° 

JJ -20.0 

2 96 +015 



300.1 



0.39+0.06 
2.13; 
0.83- 



t +0.14 
-0.20 
; +0.10 
-0.11 



3.1 .1 . Error estimates for the Orbital Elements 



The Levenberg-Marquardt algorithm also computes the covari- 
ance matrix of the fitted orbital parameters, which allows deriva- 
tion of standard error estimates for the parameters. However, the 
algorithm assumes that the^ 2 function can be approximated by a 
quadratic form in the region close to the minimum. While this is 
a reasonable assumption for iteratively searching the minimum, 
it is clearly a rather bad assumption if one is interested in confi- 
dence limits for the fitted parameters (cf. Fig.[3]and|4|. Better er- 
ror estimates c an be obtained b y studying the^ 2 function around 
its minimum dPress et alj|1992l) . Since we are interested in the 
confidence interval for each parameter taken separately, we have 
to perturb one parameter (for example To) away from the min- 
imum, and optimize all the other parameters. Any perturbation 
of a parameter will of course lead to a larger^ 2 . The range in 
To within which x 2 (To) - X^m < ^ defines the 68 % confidence 
interval for To- This interval is usually not symmetric around 
the To of the best fit, therefore we list in Table [3] separate limits 
for positive and negative perturbations. It should be noted that 
these limits describe the parameter range that contain 68 % of 
the probability distribution, which is equivalent to the commonly 
used lcr-errors. However, the errors are not normally distributed, 
therefore a 2cr interval will not contain 95 % of the probability 
distribution. 

Orbits with x 2 ~ xtm < 1 are not necessarily all within a 
single region in parameter space around the global minimum. 
There may be local minima that also fulfill the criterion. In fact, 
some of our Levenberg-Marquardt fits converged on orbital so- 
lutions with eccentricity e > 1 and x 2 - Xmi n < ^ ■ These orbits 
are within the 68 % confidence limit. Therefore, we cannot fully 
exclude that the system is unbound, although it is more likely to 
be bound. 

The system mass Ms given in Table[3]is computed from the 
semi-major axis and the period. To convert the angular separa- 
tion in milli-arcseconds (mas) to a linear separa tion in AU, a 
distanc e to T Tau of 146.7 ± 0.6 pc was adopted (Loinarc fet alj 
l2007bf l The semi-major axis a and the period P of the T Tau S 
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3 lLoinard et all {2007b) give a distance of 147.6 ± 0.6 pc in the ab- 
stract of their paper, but 146.7 ± 0.6 pc in the main text. Only the latter 
is compatible with the parallax they report. 



-50 -100 

Offset from Sa in R.A. [mas] 

Fig. 6. Best fitting orbit of T Tau Sb around Sa, in the rest frame 
of Sa. The observed positions are marked by their error ellipses 
and lines connecting the observed and calculated position at the 
time of the observations. The dash-dotted line indicates the line 
of nodes, the dashed line the periastron, and the arrow shows the 
direction of the orbit al motion. The two da sh-dotted ellipses are 
the measurements by Mavam aet alj d2006l) . 



binary are correlated around the minimum of the ^-surface: If 
one sets the period to a value higher than at the minimum^- 2 and 
optimizes the other 6 parameters, then the resulting a will also 
be larger than at the global minimum. This has to be taken into 
account in the computation of the error of the system mass. One 
way to do this is to take into account the covariance of a and P. 
However, since the covariance matrix is a good approximation 
only in a small region around the^ 2 minimum, we choose a dif- 
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Fig. 7. Detailed view of the orbit of T Tau Sb in the rest frame of 
Sa, and the observations since 2002. The symbols are the same 
as in Fig. [6] 



ferent approach: During the determination of the confidence in- 
tervals for a and P, we also compute the system mass. The max- 
imum and minimum masses within the region where A^ 2 < 1 
are adopted as confidence limits and listed in Table [3] It is worth 
noting that we did not have to decide whether we want to use the 
results from the confidence interval in a or in P - the orbits are 
identical. This is caused by the strong correlation between a and 
P, which in turn leads to a small error in Ms compared to the 
larger uncertainties in a and P (see also Sect. 14.21 ). 



3.2. The Orbit of T Tauri S and the Mass Ratio Sb/Sa 

The positions of T Tau Sb relative to Sa allow us to determine 
only the combined mass of Sa and Sb. To compute the individ- 
ual masses, we need to know the mass ratio q, which can be 
computed if the position of the center of mass (CM) of Sa and 
Sb is known. Unfortunately, we cannot observe the CM directly. 
However, we know that the CM of Sa and Sb is in orbit around 
T Tau I'fl, and that Sa and Sb are in orbit around their CM (see 
Fig.©. 

One might think that in order to derive the mass ratio q, it 
is necessary to create a model of the entire triple system, and fit 
both the orbit of the CM, the orbit of Sa and Sb, and the mass 
ratio simultaneously. However, we know that the CM is always 
on the line between Sa and Sb, and its distance from Sa is the 
constant fraction q/(l + q) of the separation of Sa and Sb. The 



4 In reality, both N and the center of mass of Sa/b are in orbit around 
the center of mass of the entire system, but this does not change the 
formalism if only relative positions are used. 
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Offset from T Tau N in R.A. [milli-arcsec] 

Fig. 8. Schematic representation of the stars and orbits in the T 
Tauri system. The origin of the coordinate system is fixed on T 
Tau N (marked by a plus sign), therefore the center of mass of 
Sa and Sb describes an orbit around N. The orbit resulting from 
our model fits is indicated by the dashed line, while the position 
of the CM of Sa and Sb on February 20, 2000 is denoted by an 
asterisk. The orbits of Sa and Sb around their CM according to 
our model fits are drawn as dotted lines, the positions of the stars 
on February 20, 2000 (computed from the model) are marked by 
a diamond and a triangle. 



positions of Sa and Sb are usually observed at the same time as 
N, so we can use the observed separation vectors Sa-Sb instead 
of a model orbit. 

Our model describes the position of the CM of Sa and Sb 
in two ways: First, it is on a Keplerian orbit around N, which 
is described by 7 orbital elements. Second, the position of the 
CM can be computed from the positions of Sa and Sb, and the 
mass ratio (which is treated as a free parameter). Standard error 
propagation is used to obtain an error estimate for this position. 
To compute^ 2 , we compare the position of the CM from the or- 
bit around N with the positions derived from the observations. 
Our model has therefore only 8 free parameters, the 7 elements 
which describe the orbit of the CM of Sa/b around N, and the pa- 
ramet er / = ql(\ + q). The parameter / is often called fractional 
mass dHeintz Il978l) . since it is the secondary star's fraction of 
the total mass in a binary. It is useful in our case because it also 
describes the fractional offset of the CM from Sa, i.e. the separa- 
tion between Sa and the CM divided by the separation between 
Sa and Sb. For a grid-search, / is better suited than q, because 
/ is confined to the range to 1 , while q is a number between 
and infinity. 

The fitting procedure is similar to that used for the orbit 
of Sa/b, except that the grid-search is carried out in 4 dimen- 
sions: eccentricity e, period P, time of periastron To, and the 
fractional mass /. Singular Value Decomposition was used to fit 
the Thiele-Innes constants, which give the remaining orbital ele- 
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merits. We stress that the orbital elements in this fit describe the 
orbit of the N/S binary, only the fractional mass / refers to the 
pair Sa/b. 

We do not expect that the elements of the orbit of S around 
N will be well constrained by the small section of the orbit ob- 
served so far. The orbit mod el is used mainly as a way to smooth 
out the measured positions. Duch ene et al.l d2006l) used second- 
order polynomials to do this, but we believe that it is better to use 
a physical model, i.e. a Kepler-orbit. With a physical model, we 
can check whether the parameters of the orbit are reasonable, 
e.g. the period, semi-major axis, and the system mass (see be- 
low). Furthermore, a Kepler-orbit is fully described by 7 param- 
et ers, only one paramete r more than the two polynomials used 
bv lDuchene etaO d2006l) . 

For the four dimensional grid search, we used 80 points for / 
in the range to 0.8, 100 points for e between and 1, 125 points 
for P in the range 100 to 3100 years, and initially 100 points for 
Tq distributed uniformly over one orbital period. Similar to the 
fit for the orbit of Sa/Sb, the grid in Tq was refined until the 
grid spacing was less than 3 days. The optimum parameter sets 
for all /, e, and P points were stored, but only the parameter 
set for the best 7q point. The most important insight we gain 
from this grid search is that the orbital parameters are not very 
well constrained, while \ 2 as a function of / shows a distinct 
minimum. 

Finally, the grid-point with the smallest^- 2 was used as start- 
ing point for a Levenberg-Marquardt fit. Since this fit method 
assumes that the measurement errors are independent of the fit- 
ting parameters, / cannot be treated like the other parameters. 
Instead, we perform a grid search over a narrow range in / to 
find the minimum. 

Fits to small sections of an orbit often result in unphysically 
high system masses (several 100 M ). Such an orbit is clearly 
not a good solution, even if its y 2 is very small. To prevent the 
Levenberg-Marquardt fit from exploring unphysical regions of 
parameter space, we added a mass constraint to the computation 
of* 2 : 



Table 4. Parameters of the best orbital solution for N - S. 



X 



z 



f i,obs ™ /,model 

Ar ;>obs 



M, 



model 



AM„, 



Here r, j0 b s and Ar l o b s are the observed position at time i and the 
corresponding error, r, mo dei is the position at time i computed 
from the model, M est and AM est are the estimated system mass 
and its error, and M mo <jei is the system mass computed from the 
orbit model. We use M est = 5 M G and AM est = 1 M G . This is 
a reasonable estimate for the total mass of T Tau N, Sa, and 
Sb, based on our r esults for Ms in Se ct. 13.11 and the mass of 
T Tau N derived bv lLoinard et al. d2007bl) . who estimated it from 
the spectral energy distribution. From comparison with theoret- 
ical pre-main-sequence tracks, they derived a mass of 1.83*2 ^9 
or 2.14^g}yM (depending on the theoretical tracks used). 

An important question is whether we want to include ob- 
servations that did not resolve T Tauri S in the fit. These mea- 
surements yield only the position of the center of light of Sa/b, 
not the center of mass. The offset between the two centers is 
unknown and not eve n constant, since b oth components are pho- 
tometrically variable dBeck et alj |2004). On the other hand, the 
data obtained before the binary nature of T Tauri S was discov- 
ered almost double the fraction of the orbit that has been ob- 
served so far. 

A fit only to observations that did resolve T Tau S results in 
an orbit model that does not reproduce the positions measured 



Orbital Element 


Value 


Date of periastron To 


2452981 tjj™ 
(2003 Dec 7) 


Period P (years) 


73940 +1 1000 
iJi ™ -12000 


Semi-major axis a (mas) 


QOQ'7 + 1 8000 

yoo i 5200 


Semi-major axis a (AU) 


1443 


Eccentricity e 


93 +0 05 

v yJ -0.05 


Argument of periastron to (°) 


198 +1 


P. A. of ascending node £1 (°) 


351 -u 


Inclination i (°) 


43 + _t 


System mass Ms (M e ) 


5 57 +o u 

-0.29 


reduced x 1 


3.3 



before 1995 (Fig.|9^). Although we cannot derive the exact po- 
sition of the center of mass, it is unlikely that the observed center 
of light is as far off as this model suggests. Furthermore, the mass 
ratio this fit yields is 2.42, i.e. according to the model, T Tau Sb 
would be 2.42 times more massive than T Tau Sa. With a system 
mass of 3 M Q (see section [370 , the mass of T Tau Sb would be 
2.1 M , inconsistent with its spectral type of about Ml. 

A compromise between using all observations and only ob- 
servations that resolved T Tau S is to use all observations, but 
add some additional uncertainty to the errors of the unresolved 
observations, in order to compensate for the offset between cen- 
ter of light and center of mass. Since the separation between T 
Tau Sa and Sb was about 50 mas in 1997, 25 mas is a reason- 
able upper limit for the separation between the two centers. In 
principle, we could use the predicted position angle of the Sa- 
Sb vector to constrain the direction of the additional uncertainty. 
However, the predicted position angle depends strongly on the 
period of the orbit of Sa and Sb around each other, which is 
rather uncertain. Therefore, we decided to not apply any con- 
straint on the position angle of the vector from center of light 
to center of mass, we simply add the same uncertainty in both x 
and y. 

Figure|9j) shows the result if 25 mas are added in both x- and 
y-direction. The resulting orbit model is closer to the unresolved 
observations than the model of Fig. [9^, within less than 2<r of 
the enlarged uncertainties. The mass ratio is 0.77, so T Tau Sb 
should have a mass of about 1 .3 M Q , which is still inconsistent 
with spectral type Ml. 

Figure [9J; shows the result of a fit where we added 10 mas 
additional uncertainty to the unresolved observations. The mass 
ratio of the best-fitting model is about 0.49, in which case the 
mass of T Tau Sb is 1 M Q . 

Finally, Fig. [9jl shows the result if we do not add any ad- 
ditional uncertainties, i.e. take all the measurements and their 
errors as they were published. The mass ratio is 0.39 + 0.06. The 
uncertainty of the mass ratio was estimated by finding the mini- 
mum and maximum q that yields ay 2 less than^^ + 1 (Fig.|T0T>. 
The lcr-uncertainty for q is half the width of this interval. With 
this mass ratio, T Tau Sb would have a mass of 0.83^' J jM Q , 
which would be in reasonable agreement with its spectral type. 
Table [4] lists the orbital elements of the best fit, together with 
estimates for their confidence limits. These confidence limits 
give the size of the region around the global minimum where 
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Offset from T Tau N in R.A. [mas] 



Offset from TTau N in R.A. [mas] 



Offset from T Tau N in R.A. [mas] 



Fig. 9. Results of fits for the orbit of the center of mass of T Tau S around T Tau N. The four panels differ in the way observations 
were treated that did not resolve T Tau S into Sa and Sb. These measurements are marked by crosses, measurements that did resolve 
T Tau S are marked by asterisks. Note that the crosses show the center of light (which can be observed directly), while the asterisks 
show the center of mass (which depends on the observations and the model-parameter q). Panel a shows the result if unresolved 
observations are not used at all in the fitting procedure. The unresolved measurements are plotted for comparison only. The orbit 
in panel b was obtained by adding 25 mas uncertainty to the unresolved observations, while panel c shows the result if 10 mas are 
added. For Panel d, the uncertainties given in the literature were used, without any additional errors to compensate for the offset 
between center of mass and center of light. In all panels, the dotted line shows the path of the center of mass of T Tau S in our 
models. 



y 2 < y 2 ■ +1. We did not search for local minima elsewhere 

^ A mm 

in parameter space. It is therefore possible that there are orbits 
with, e.g., significantly different inclination, which also fulfil the 
X 2 < + 1 criterion. In particular, x 2 as function of period 
or semi-major axis is rather flat and appears to be dominated 
by random fluctuations, making it difficult to find a meaningful 
confidence limit. The orbital elements given in Table [4] should 
be good enough to predict the positions up to a few decades in 
the past or future (depending on the required accuracy). 

We conclude that the mass ratio Ms^/M^ is probably in the 
range 0.3 - 0.5, but the time-span of observations that resolve 
the Sa-Sb-pair is not long enough to exclude higher mass ratios. 

Figure QT] shows the motion of T Tau Sa and Sb in the rest 
frame of T Tau N, including both the observed positions and 
those predicted by our fits. 



4. Discussion 

4.1. How realistic is our model? 

Our orbit model assumes that the center of mass of T Tauri S is 
on a true Keplerian orbit, which is not disturbed by the binary 
nature of the southern component. Whether this assumption is 
valid is difficult to judge with the data at hand, since it depends 
on the minimum distance of Sa and Sb to N, which in turn de- 
pends on the shape and orientation of both orbits. 

Our best fit for the orbit of T Tau Sa and Sb around each 
other has a maximum separation of about 46 AU. The projected 
separation between T Tau N and S is currently about 100 AU. 
This would be close enough to the separation of Sa and Sb to 
cause significant perturbations on the orbit of Sa/Sb. However, 
we know only the projected separation, little is known about the 
3-dimensional distance. The best orbit resulting from the fits de- 
scribed in Sect. 13.21 has a semi-major axis of about 1500 AU, 
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Fig. 11. Motion of T Tau Sa and Sb in the reference frame of 
T Tau N. The path of Sa predicted by our orbit models is shown 
by the dotted line, the path of Sb by the dashed line, and the 
path of their center of mass (CM) by the dash-dotted line. The 
predicted positions on January 1, 1997 to 2008 are marked by 
crosses. The observed positions of Sa are depicted by diamonds, 
those of Sb by triangles. The solid lines connect the observed 
and the predicted positions. 



albeit with a high eccentricity, which leads to a periastron dis- 
tance of 105 AU. In fact, the semi-major axes of N-S orbits with 
a small y 2 s P an the range from about 180AU to 1800AU, but 
all their periastron distances are between 100 and 110AU. This 
might indicate that the orbit of Sa and Sb around each other is 
indeed perturbed by T Tau N, but we have to wait for obser- 
vations of a much larger fraction of the N-S orbit to draw any 
conclusions. With the data currently available, it is not possible 
to detect deviations from a Kepler-orbit. 

4.2. Why is the mass so well constrained? 

It is surprising that the system mass of T Tauri S derived from 
our orbital fits has a much smaller fractional error than the semi- 



major axis and the period. Mathematically, this follows from the 
correlation of a and P around the minimum^ 2 . This correlation 
is found as a result of the fitting procedure. However, there is 
also a physical reason for it. 

Newton's law of gravitation states that the acceleration vec- 
tor r of the companion is directed towards the primary star, and 
that it is proportional to the system mass M and the inverse of 
their distance r squared: 

r = GMr 2 , 

where G is the constant of gravitation. If we could measure r 
and r at any point of the orbit, we could compute M. However, 
we cannot observe the component of r along the line of sight. 
There are only two points in the orbit where this component is 
zero: the nodes of the orbit. At these points, we can observe r 
and r and compute the mass. In fact, according to our best or- 
bit model, T Tauri South passed through a node in 2002, in the 
middle of the time interval covered by the observations. In prin- 
ciple, we can use the observations near this time to compute the 
mass directly from r and r. In practice, uncertainties in the mea- 
surements prevent us from obtaining a reliable result. However, 
any orbit model approximating the measured positions will also 
contain good approximations for r and r, and therefore lead to 
the same system mass. This means that the system mass does not 
vary much among orbit models approximating the observations. 
A consequence of this nearly constant mass is the correlation of 
a and P, since Kepler's third law states that a 3 is proportional to 
P 1 , if the mass is constant. 

4.3. Comparison to previous results 

IDuchene et al.l (120061) presented an orbital solution for the T Tau 
S binary based on IR measurements published before 2006 and 
radio observations at 2 and 3.7 cm wavelength (they assume that 
the radio source is identical to T Tau Sb). The y adopted the 
distan ce to T Tau of 141.5 ± 2.8 pc estimated by lLoinard et al.l 
(2005), while we use t he more recent value of 147.6 ± 0.6 pc 
(lLoinard et alJ l2007bh . IDuchene et alJ (120061) find a period of 
21.66 + 0.93 yrs, a semi-major axis of 82.1 ± 1.8 mas, and a 
total mass of T Tau S of 3.79 + 0.41 M (scaled to the more 
recent distance determination). This is inconsistent with our re- 
sults, and they estimate uncertainties that are almost two orders 
of magnitude smaller than ours. We believe that this is caused by 
the method they used to d erive the uncertainties. 

IDuchene et al.l d2006h employed a bootstrap method to esti- 
mate the errors in the fitted orbital elements. They created artifi- 
cial datasets where the data points were randomly drawn from 
Gaussian distributions centered on the actual measurements. 
The ^-minimization routine was run on each of the artificial 
datasets, which results in a distribution for each of the parame- 
ters. The widths of these distributions were used as an estimate 
for the errors of the fit parameters. 

This bootstrap method is not a good way to estimate errors 
in the case of T Tauri S. We have only observational data for 
about one quarter of the orbit. Therefore, a wide range of, e.g., 
periods is compatible with the data (see Fig. [3]). In the part of 
the orbit that has been observed so far, all these orbits are close 
to the observed data points, and hence also close to each other. 
The orbits show very different separations in the parts that have 
not (yet) been observed. The uncertainty in the orbital param- 
eters is caused by incomplete coverage of the orbit, not by the 
measurement errors of the observations. Bootstrapping methods 
help in determining how errors in the input data propagate to 
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errors in the results. In the case of T Tauri, the measurement er- 
rors contribute only a small part of the uncertainty of the orbital 
elements. Therefore, bootstrapping methods underestimate the 
errors. 

The best way to estimate confidence limits on the fitted pa- 
rameters in this case is to ana lyze y 2 as a function of the param- 
eter values dPress et al.l lT992, see also Sect. 13. It . The best esti- 
mate for the fitted parameters is, of course, the point where x 2 is 
a minimum. As confidence limit for each parameter, we use the 
interval within which^f 2 -^ 2 nin < 1 ■ These limits contain 68 % of 
the probability distribution. A note of caution is indicated here: 
The orbital elements are most certainly not normally distributed. 
Therefore, a single <x is not sufficient to describe the distribu- 
tion. First, the confidence limits are not symmetric around the 
minimum (cf. Table [3). Second, the interval containing 95 % of 
the probability distribution is not twice as large as the interval 
containing 68 %. The well-known rules for normal distributions 
do not apply here. The 95 % confidence interval is defined by 
X 1 - Xmin < 4, but since x 2 is not a quadratic function of the or- 
bital parameters, the interval is not twice as large as the interval 
defined by x 2 ~ xtnm < 1 • We determined only the 68 % confi- 
dence intervals and list them in Table [3] they should suffice to 
indicate the uncertainties in the parameter estimates. 

5. Summary 

We have collected all available near-infrared astrometric data on 
the T Tauri system from the li terature. We also pre sent a new 
analysis of the data published in Kohle ret al.l 12000a) and 5 new 
data points obtained with NACO at the VLT. Binary orbit models 
were fitted to the relative positions of T Tau Sa and Sb in order to 
estimate the orbital elements. We find that most elements are not 
very well constrained, e.g. the period is 93*il years. The total 
mass of the T Tau S binary, however, can be estimated more pre- 
cisely to be 3.0^o 24 M Q . It is unlikely that T Tau Sb was recently 
ejected from the system and is now in a highly eccentric orbit or 
even escaping from the system. However, some orbit solutions 
with e > 1 are within 1 cr of the best fit, so an unbound orbit 
cannot be ruled out completely. 

We used the positions of T Tau Sa and Sb relative to T Tau 
N to fit a combined model of the N-S and the Sa-Sb binary. The 
separation of Sa and Sb is taken from the observational data, 
only the mass ratio Msb/Msa is a free parameter for the fit. Due 
to the small fraction of the N-S orbit observed so far, no reliable 
constraints for the orbital elements of the N-S binary could be 
obtained from the data, but the mass ratio of the southern binary 
can be estimated to be 0.4 + 0.1. This corresponds to individual 
mass estimates for T Tau Sa and Sb of 2.1 +0.2 and 0.8+0.1 M B , 
The spectral type of T Tau Sb was used to decide which orbital 
solution is acceptable, it is therefore not surprising that the mass 
estimate is in agreement with the spectral type. 

Our results indicate that T Tau Sa is at least as massive as T 
Tau N, although it is much fainter at optical wavelengths. This 
can be explained by extinction due to circumstellar (or circumbi- 
nary) material, possibly an edge-on disk. Observations at longer 
wavelengths with even higher spatial resolution are necessary to 
study this material (Ratzka et al., in prep.). 
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